Nucleation in finite topological systems during continuous 
metastable quantum phase transitions 
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Finite topological quantum systems can undergo continuous metastable quantum phase transitions 
to change their topological nature. Here we show how to nucleate the transition between ring 
currents and dark soliton states in a toroidally trapped Bose-Einstein condensate. An adiabatic 
passage to wind and unwind its phase is achieved by explicit global breaking of the rotational 
symmetry. This could be realized with current experimental technology. 



Phase transitions have long been considered equilib- 
rium phenomena of infinite systems that can involve in- 
teresting non-equilibrium nucleation dynamics, e.g. the 
formation of topological defects [I]. In recent years, 
quantum phase transitions have also been identified in 
excited states of nuclei as well as ultracold quantum gases 
[2H5] by the nonanalytic change in spectral properties, or 
their finite-system precursors, as a system parameter is 
changed. Such phase transitions are not manifested in 
equilibrium but they have dynamical consequences 0] . 
So far, little is known about the nucleation process. 

Nucleation of equilibrium quantum phase transitions 
often relies on external symmetry breaking, which can 
occur locally. An example is the Ising model [6], where a 
tiny external magnetic field is sufficient for breaking the 
rotational symmetry of the spins and allowing the system 
to equilibrate. Similar mechanism were discussed in the 
context of vortex nucleation in Bose gases [7HH] ■ Here 
we theoretically investigate the quantum phase transi- 
tion between metastable states of a Bose gas in a ro- 
tating toroidal trap identified in Ref. [5]. The transi- 
tion occurs between topological vortex states and non- 
topological soliton states. Vortex states of N atoms cor- 
respond to ring currents and carry integer angular mo- 
mentum per particle L/(Nh) = J. As the rotation fre- 
quency of the trap fl is changed, the ground state jumps 
between vortex states [TU] and thus the integer J is a 
topological charge. The soliton states carry non-integer 
L/(Nh). In mean-field theory they are approximated by 
dark solitons [11], which carry a localized density notch 
and thus break rotational symmetry in addition to chang- 
ing the topological nature of the system, although no 
symmetry breaking is required for the phase transition 
in the quantum description of the finite system. 

In this Letter, we show how to nucleate the phase tran- 
sition by means of a global symmetry breaking potential, 
which creates an adiabatic passage through metastable 
states. Although symmetry can be restored for a finite 
system, the passage naturally leads to the emergence of 
symmetry broken soliton solutions for large particle num- 
ber N. In order to maintain metastability, we further find 



that local symmetry breaking needs to be avoided as it 
would lead to the unwanted thermalization of the ex- 
cited states [IS]. This is in stark contrast to equilibrium 
phase transitions, where thermalization is desirable. The 
explicit global symmetry breaking is achieved by tilting 
the trap axis by some angle 9 and rotating it with the 
frequency fi. The frequency f2 is decreased from a max- 
imum value ^initial to a minimum value f2fi na i within a 
finite time interval, after which the tilt is decreased. Un- 
der certain conditions angular momentum is established 
in the Bose gas up to a dark soliton state or a vortex 
state as seen in Figs. [T] and [2] respectively. The protocol 
is counterintuitive as we decrease the rotation frequency 
of the system to increase its angular momentum. It is 
reminiscent of the stimulated Raman adiabatic passage 
used in quantum optics to populate a metastable state 
through an intermediate level by a counterintuitive pulse 
sequence [T3] . 

Interacting bosonic atoms of mass m confined to a ro- 
tating toroidal trap are described by the Hamiltonian 
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where ^> = ^(x) is the atomic field operator, U the 
bare interaction constant assumed repulsive (U > 0), 
and e = mgRsm(6) the strength of the symmetry break- 
ing potential (where mg is the gravitational force). The 
toroidal trap is approximated by a one-dimensional ring 
of the radius R. 

Let us first discuss the outlined problem in the mean- 
field approximation. Using a rotating coordinate frame 
with frequency £1 and introducing the scaled variables 
(j> = x/R, t = tEo/h, lo = nil/(2E ), the Gross- 
Pitaevskii (GP) equation reads [IT] 
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where 7 = U(N - l)/(2nR) and E = H 2 /(2mR 2 ) 
gives the relevant energy scale. The classical field x 



is normalized as L 27r |x| 2 d</> = 1 and it is periodic on 
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<p, x(0) = x(27r). We apply Broyden's method [2] to 
solve Eq. ([2]) numerically for the stationary solutions 
x(t, 4>) — x(^) e_VT ) where /i-Bo is the chemical poten- 
tial. For e = 0we present in Fig.^a) three solutions cor- 
responding to two vortex states x(4>) — e lJ */v / 27r with 
J = and J = 1, which are connected by a soliton 
branch. Finite e opens a gap in the soliton branch, see 
Fig. ^b) . The lower branch corresponds to a soliton sit- 
ting on the crest, while the upper branch connects to a 
soliton located in the trough of the external potential. 
In the following we will show that the lower branch is 
dynamically unstable, while the upper branch is stable. 
Thus an adiabatic passage through the upper branch is 
possible. Numerical simulations of the system dynamics 
on the basis of Eq. ^ confirm this expectation. The 
solid lines in Figs. |ljb) and [2jb) show the momentum 
per atom as a function of time for the adiabatic proto- 
cols shown in the top panels. The final soliton and vortex 
states are shown in Figs. [2jc) and (d), respectively. 

Mean-field theory thus presents the following picture: 
The bifurcation points between the soliton and vortex 
branches seen in Fig. [3^ a) mark the transition between 
rotationally symmetric, and symmetry broken [Fig. [2jc)] 
phases. With finite e, the symmetry is formally broken 
everywhere and the bifurcation makes way to a swal- 
lowtail structure [TS] in Fig. [3jb) . In contrast to the 
dynamics of continuous phase transitions in infinite sys- 
tems where adiabaticity is always violated pp, an adia- 
batic passage is established (arrows in Fig. [3| to change 
the symmetry properties of the initial state. Due to the 
global nature of the symmetry broken state with a single 
soliton in the whole ring there is no formation of domain 
structures with locally broken symmetry in further con- 
trast to Refs. pp. 

We proceed with discussing and quantifying the nec- 
essary conditions for an adiabatic passage. Useful in- 
sight into the physics of the considered process is ob- 
tained within a two-mode approximation, which is jus- 
tified for small 7 <C Eq. Neglecting all coefficients in 
the Fourier expansion except k — 0,1 of the function 
X, x{ T i4>) — Sfc^fc( T ) exp(i/c</>), the system dynamics is 
described by the effective Hamiltonian 
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where I = |6i| 2 and 1? is the relative phase for the am- 
plitudes 60 and b\. The phase portrait of ^ contains a 
stability island around the elliptic point (/, 1?) = (/*,0) 
(cf. Fig.flj:, d), where I* « 1/2 + (1/2 + lo)E /j depends 
linearly on u) and reaches values of zero and unity for 
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respectively. We note that in the lab frame this stability 
island corresponds to a nonlinear resonance. Thus the 
adiabatic passage has a simple physical interpretation: 




FIG. 1. Adiabatic passage from the ground to the dark soli- 
ton state, (a) Protocol for variations of the perturbation pa- 
rameter e (dashed line) and the precession frequency uj (full 
line) for 7 = 0.8i?o. (b) Angular momentum per particle 
in the mean-field approximation ([2| over time (full line) and 
probability density of the quantum two-mode simulation with 
Hamiltonian Q on the same axes. White and black corre- 
spond to zero and maximum probability density, respectively. 
Phase portraits of Eq. |3| are presented for ui — 0.7 (c) and 
u = 0.5(d). 



by ramping e to a finite value we capture the system 
into the nonlinear resonance, transport it to any desired 
value of / by adiabatically changing the frequency u from 
^initial ^+ to Wfi na i, and release the system by ramping 
e back to zero. The necessary condition for the time scale 
of this process is At 3> fij 1 , where O s « y/je/h is the 
frequency of small oscillations near the elliptic point of 
the stability island (the approximate sign is replaced by 
an equal sign for uj — 1/2). 

For stronger nonlinearity, 7 > Eq, when the two-mode 
approximation is not justified, the adiabaticity conditions 
and stability can be studied using the Bogoliubov ap- 
proach, which linearizes the time dependent GP equation 
([2]) [TT]. This leads to the eigenvalue problem 

XiUi = {H G p - M + ^1 / E o\x?) u i + ^l/E x 2 Vi, 



-XiVi = {H G p 
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Here, x = X^(<^) is the stationary solution for the up- 
per soliton branch in Fig. [3^b) and Hqp is given in the 
square brackets in Eq. ([2]). For the lower branch we find 
a single imaginary eigenvalue A, which indicates that the 
solutions on this branch are dynamically unstable [llj . 
The Bogoliubov analysis for the upper branch indicates 
that it is dynamically stable and that there is a gap to 
the lowest lying elementary excitation. The correspond- 
ing frequency fi s = XE /h has the asymptotic behavior 
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FIG. 2. Adiabatic passage from the ground state to the J = 1 
vortex state, (a) Protocol for variations of the perturbation 
parameter e (dashed line) and the precession frequency uj (full 
line) for 7 = 0.8i?o. (b) Angular momentum per particle 
in the mean-field approximation (full line) and probability 
density of the corresponding quantum simulation as in Fig. [I] 
The dark soliton state (c) and the vortex state (d) represent 
the final states of the corresponding mean field simulations. 
Full lines represent density, dashed lines represent phases. 



For small nonlinearities 7 <C Eq this result agrees with 
the previous two-mode analysis. 

The large 7 result can be understood by noting that the 
soliton is a localized object. In this case fl s is interpreted 
as the frequency of small oscillations of the soliton around 
a stationary point. Indeed, it can be shown that the small 
oscillations of the soliton around the stationary points are 
described by the equation |16j 



dt 



± 



2i? 2 



const , 



(7) 



which is valid for ft w Eo/h. The ± sign corresponds 
to a soliton sitting at the trough or the crest of the po- 
tential, respectively. While the minus sign describes an 
unstable situation (lower soliton branch in Fig. ^p), the 
plus sign describes harmonic oscillations of a particle with 
mass 2m. The soliton located in the trough is thus dy- 
namically stable, with the frequency of small oscillations 
n s = ^e/(2mR 2 ) = VWe/h- 

In practice, the adiabaticity condition At ^> fi" 1 dis- 
cussed above should be made even stronger by requiring 
that accelerated rotation should displace the soliton by 
less than the the soliton core size, which is of the order 



of the healing length £ = 
condition 



h 2 /2mj. This yields the 
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FIG. 3. The energy of the GP equation (upper row) and 
the energy spectrum of the Hamiltonian (|9| in the 2-modes 
approximation (lower row). Parameters are 7 — O.8E0, and 
e = (left column) and e = 0.04_Eo (right column). The 
arrows show the passage under adiabatic change of the fre- 
quency u>. Three distinct solutions exist in panel (a) between 
the bifurcation points u- w 0.13 and lj+ w 0.87. 



firms the estimate ([8]) . For the parameters of Fig. [T] the 
right hand side of ^ifis « 0.04(i?o/?i) 2 . It is seen in Fig.[Tj 
and[i]that for a ten times smaller rate (l = 0.004(I?o/7i) 
the dark soliton and the vortex states (green solid lines) 
are reached. 

The mean-held analysis presented so far can be sub- 
stantiated by a quantum analysis. Using scaled variables, 
a rotating coordinate frame, and expanding the held op- 
erators in the angular-momentum basis, the Hamiltonian 
of the system takes the form of the Lieb-Liniger Hamilto- 
nian |17j with an external symmetry-breaking potential 
under periodic boundary conditions 

e 
2 
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We are interested in the time evolution with initial condi- 
tion given by the nonrotating ground state (vortex state 
with J = 0). The soliton states have dehnite angular 
momentum for e — and, within the validity of the two- 
mode approximation, are given by 



\in) = (bl) N - l (b\) l \y a c) 



(10) 



The numerical simulation of the system dynamics con- 



Truncating the Fock basis to the states ([10b, the Hamil- 
tonian (|9| is a tri-diagonal (N + 1) x (N +1) matrix. 
Fig. [3^c) shows the spectrum of this matrix for N = 10, 
7 = O.8S0, e = 0, and < u < 1. Note that for U = all 
levels would cross at one point at uj — 1/2. Finite interac- 
tions remove this degeneracy, leading to the appearance 
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of a caustic in the level crossing pattern. The spectrum 
of the system for a finite e = 0.04i5o is shown in Fig.[3]jd) 
[18] . Now all level crossings are substituted by avoided 
crossings. The Hamiltonian in the two mode model can 
be recast to take the form of a quantized version of the 
the effective Hamiltonian ([3|, where the action variable 
I is now associated with the operator / = —i(l/N)d/di9 
and 1/N plays the role of the Planck constant. In partic- 
ular, we find that the transition frequency between the 
caustic levels is given by f2 s of Eq. Q . 

Within the quantum model we can find a condition for 
the frequency Wi n iti a i from the crossing point of the two 
energy levels that are lowest for w = in Fig.[3|c). In the 
two mode approximation, we obtain lu + of Eq. Q, the 
same result as from the mean-field analysis. A more thor- 
ough comparison of the upper and lower rows in Fig. [3] 
indicates that energies of the soliton states calculated 
within the mean-field approximation are shifted in the 
negative direction compared to the quantum calculation. 
This is a manifestation of the occupation of additional 
angular momentum modes outside the two-mode model. 
In order to extend the result for w+ beyond the two- 
mode approximation, we consider the rotation frequency 
uji = y/l 2 + 2~//E /2 where the Bogoliubov phonon with 
angular momentum lh > acquires zero energy. The ex- 
pression for uj + = lu\ thereby generalizes Eq. Q to larger 
values of 7. Generally, it looks possible to excite the vor- 
tex state with topological charge J with the condition 
luj < ^initial < |23j through metastable states with 

J dark solitons [5]. 

It is interesting to simulate the dynamics of the quan- 
tum system ([9| for the same protocols as were used in 
the mean-field simulations. Due to exponential prolif- 
eration of the Hilbert space with N, this can be done 
only for a small number of atoms N ~ 10. The density 
plots of Figs. |TJb) and [2jb) show the dynamics of the 
system Q in the two-mode approximation for N = 10 
and 7 = 0.8i?o- (We have checked that for these value 
of the interaction constant and number of particles the 
result remains unchanged if we use a four-mode approxi- 
mation.) The gray scale encodes populations of the Fock 
states (JToJ) . It is seen in Fig. [I] and [2] that our protocols 
almost entirely populate the target excited state \tpN/2) 
(Fig. [jjb)) and \ip N ) (Fig. ^b)), respectively with the 
probabilities 0.93 and 1. 

The quantum model ([£]) reveals why global symme- 
try breaking is effective for nucleating the metastable 
phase transition: The long wavelength perturbing poten- 
tial proportional to e couples neighboring soliton states 



(10 1 effectively and thus creates the avoided crossings 



that provide the adiabatic passage. A short wavelength, 
local symmetry breaking perturbation, in contrast, would 
tentatively couple to excited states with larger angular 
momentum difference, generate drag, and thus compro- 
mise the adiabatic passage. For a large system (R ^> £) 
this is expected from both quantum [10] and mean-field 



[TU] consideration of a localised impurity moving at su- 
personic velocity ftR > fliR = c + 0(R^ 1 ), where 
c = \f~fjm is the speed of sound. 

There is still a crucial difference between the quan- 
tum and mean field models: While the soliton states of 
the GP approximation break the rotational symmetry, 



the corresponding quantum states (10) do not. Instead, 



they correspond to fragmented condensates, where both 
the occupation of the k — and the k = 1 mode can 
become large. Whether we obtain a fragmented conden- 
sate with preserved rotational symmetry or a condensate 
with broken symmetry, depends on the rate for changing 
the parameter e at the very end of the passage. Simple 
calculation reveals that adiabatic condition is fulfilled if 
the rate satisfies de/dt < ~/ 2 /(2hN) 20J. For large N 
it becomes more difficult to restore symmetry, since the 
time scales required for restoring the symmetry become 
large. When e is changed back to zero in finite time, the 
system 'sits' on the caustic where the density of states 
tends to infinity if N — > 00. That a large particle num- 
ber brings a qualitative change is familiar from the gen- 
eral principle of symmetry breaking in condensed matter 
[21) . In contrast to conventional phase transitions, how- 
ever, where symmetry breaking occurs spontaneously, the 
metastable quantum phase transition discussed in this 
Letter requires an explicit breaking of its symmetry. 

So far we have mostly considered the case of moder- 
ate nonlinearity, where the two-mode approximation is 
justified. In typical experiments with BECs, larger non- 
linearities beyond the validity of the two-mode approxi- 
mation are relevant. We have verified within the mean- 
field theory that adiabatic passage for stronger nonlin- 
earities is possible as well [23] ■ We estimate the time 
scale for the adiabatic passage to a dark soliton state us- 
ing Eq. Q and setting e ~ 7 as At > mR 2 /h. In the 
recent experiment with 23 Na [22 , R ~ 20fim. This gives 
At 3> Is. For lighter atoms, e.g. 7 Li, and smaller radius, 
e.g. R ~ 10/zm, one can achieve the condition At 3> 0.1s 
and 9 « 0.5°. Therefore, time scales of the order of Is 
will provide an adiabatic passage. 

In conclusion, we have proposed an experimental 
scheme to nucleate a metastable quantum phase tran- 
sition in an ultra-cold Bose gas by an adiabatic pas- 
sage. We have shown that this is achieved by an explicit 
breaking of the global rotational symmetry. The pro- 
cedure generalizes nucleation procedures, which change 
the ground state symmetries in the course of ordinary 
continuous phase transitions. 
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